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[i]  Energetic  very  low  frequency  (VLF;  frequencies  <0.004  Hz)  surf  zone  eddies  (SZEs) 
were  observed  on  a  beach  composed  of  shore-connected  shoals  with  quasi-periodic 
(~125  m)  incised  rip  channels  at  Sand  City,  Monterey  Bay,  California.  Incident  waves 
consisted  of  predominantly  shore -normal  narrow-banded  swell  waves.  SZEs  were  located 
outside  the  gravity  region  in  alongshore  wave  number,  ky ,  spaced  within  the  VLF  band, 
and  did  not  appear  to  exist  in  higher- frequency  bands.  The  SZEs  were  significant 
(Un ns, vlf  ~  0.25  m/s)  and  constant  in  intensity  within  the  surf  zone  (shore-connected 
shoals  and  rip  channels)  and  rapidly  decreased  offshore.  The  alongshore  and  cross-shore 
SZE  velocity  variances  were  similar  in  magnitude.  VLF  SZE  velocities  were  not  forced  by 
VLF  surface  elevations  and  were  not  well  correlated  with  rip  current  flows  (r2  =  0.18). 

There  is  an  indication  that  the  SZEs  were  related  to  wave  forcing,  with  the  SZEs 
statistically  correlated  with  incoming  sea-swell  wave  height  (r2  =  0.49).  F-ky  spectral 
estimates  illustrate  a  strong  relationship  between  rip  channel  spacing  and  SZE  cross-shore 
velocities  (ky  =  ±0.008  m_1)  and  minimal  SZE  alongshore  velocity  variation  {ky  =  0  m-1). 

Data  analysis  suggests  that  the  SZEs  are  not  simply  instabilities  of  an  unstable  rip  current 
jet.  A  simple  conceptual  model  suggests  that  SZE  f-ky  spectra  can  be  explained  by  the 
entire  rip  current  circulation  cells  oscillating  predominantly  in  the  cross  shore  and  slightly 
in  the  alongshore.  INDEX  TERMS:  4546  Oceanography:  Physical:  Nearshore  processes;  3020  Marine 
Geology  and  Geophysics:  Littoral  processes;  4512  Oceanography:  Physical:  Currents;  4599  Oceanography: 

Physical:  General  or  miscellaneous;  KEYWORDS:  rip  currents,  surf  zone,  nearshore,  eddies,  circulation,  waves 
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I.  Introduction 

[2]  Velocity  measurements  within  the  surf  zone  for 
field  experiments  on  planar  and  barred  beaches  illustrate 
that  low-frequency  (infragravity)  motions  are  ubiquitous. 
Infragravity  (0.004-0.04  Hz)  motions  are  often  the  dom¬ 
inant  energy  on  the  beach  compared  to  the  energy  within 
the  sea-swell  frequency  band  [ Suhayda ,  1974;  Huntley , 
1976;  Huntley  et  al.,  1981;  Holman ,  1981;  Wright  et  al. , 
1982;  Guza  and  Thornton ,  1985],  even  for  a  beach  with 
persistent  rip  channels  [A.  J.  H.  M.  Reniers  et  al., 
Modeling  infragravity  motions  on  a  rip-channel  beach, 
submitted  to  Coastal  Engineering ,  2004,  hereinafter 
referred  to  as  Reniers  et  al.,  submitted  manuscript,  2004; 
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connected  shoal  beach,  submitted  to  Marine  Geology , 
2004,  hereinafter  referred  to  as  MacMahan  et  al.,  submit¬ 
ted  manuscript,  2004].  Infragravity  energy  is  confined  to 
a  region  within  f-ky  space  defined  by  the  zero  mode  edge 
wave  (/  >  f0\  \ky\  <  kyo)  (Figure  1). 

[3]  Energetic  vorticity  motions  (nongravity  motions), 
typically  occurring  at  frequencies  less  than  infragravity 
motions,  coexist  in  the  surf  zone.  The  distinction  between 
gravity  waves  and  vorticity  motions  is  determined 
through  frequency-wave  number  {f-ky)  spectral  estimates. 
We  define  nongravity  motions  for  frequencies  less  than 
0.004  Hz  as  very  low  frequency  motions  (VLFs)  since 
infragravity  waves  in  the  literature  have  traditionally  been 
referred  to  as  low-frequency  motions.  Furthermore,  VLFs 
are  considered  to  be  predominantly  composed  of  surf 
zone  eddies  (SZEs),  owing  to  the  relatively  small  space 
within  the  gravity  restoring  region  for  frequencies  less 
than  0.004  Hz.  (Figure  1). 
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Figure  1.  F-ky  spectrum  illustrating  gravity  and  nongravity  region  separated  by  mode  0  edge  wave.  The 
very  low  frequency  (VLF)  motions  are  further  defined  in  the  nongravity  region  for  frequencies  less  than 
0.004  Hz. 


[4]  Most  nearshore  studies  have  associated  SZEs  with 
shear  instabilities  of  the  longshore  current.  In  1989,  Oltman- 
Shay  et  al  [1989]  discovered  energetic  motions  in  the  long¬ 
shore  current  located  within  the  nongravity  regions  in  f-ky 
space  (Figure  1),  ranging  in  frequencies  from  0  to  0.045  Hz. 
These  motions  were  most  energetic  for  frequencies  less  than 
0.004  Hz  (i.e.,  VLFs),  decreasing  in  energy  with  higher 
frequencies.  The  longshore  current  instabilities  have  a  near- 
linear  dispersion  relationship  in  f-ky  space  related  to  the  shear, 
speed,  and  direction  of  the  longshore  current  [ Bowen  and 
Holman ,  1989;  Oltman-Shay  et  al ,  1989;  Noyes  etal .,  2004]. 

[5]  Field  investigations  of  topographically  controlled  rip 
currents  [Shepard  and  Inman,  1950;  Sonu ,  1972;  Aagaard  et 
al .,  1997;  Smith  and  Largier ,  1995;  Brander ,  1999;  Brander 
and  Short ,  2000,  2001 ;  MacMahan  et  al ,  2004]  documented 
rip  current  velocity  pulsations  within  the  infragravity  fre¬ 
quency  band.  In  addition  to  the  energetic  rip  current  pulsa¬ 
tions,  energetic  VLF  motions  have  also  been  observed 
within  rip  current  systems.  Smith  and  Largier  [1995]  acous¬ 
tically  measured  rip  current  motions  outside  the  breaker  zone 
and  found  that  they  fluctuated  aperiodically  with  temporal 
scales  greater  than  15  min.  Haller  and  Dalrymple  [2001] 
documented  oscillations  of  the  rip  current  jet-like  motions  in 
a  laboratory  study,  referred  to  as  rip  current  instabilities. 

[6]  Most  previous  field  experiments  have  either  over¬ 
looked  the  amount  of  energy  associated  with  VLFs  by  not 
sampling  sufficiently  long  to  resolve  them  or  attributed 
these  signals  to  infragravity  motions.  Brander  and  Short 
[2001]  reexamined  their  original  observations  by  Brander 
[1999]  and  found  VLFs  in  their  measurements.  Ogston  and 
Sternberg  [2003]  re-analyzed  measurements  obtained  by 


Beach  and  Sternberg  [1988]  in  the  surf  zone,  using  longer 
records,  and  found  that  there  was  a  significant  amount  of 
energy  in  the  VLF  band  for  both  velocity  and  sediment 
concentration  measurements.  Thus  VLFs  also  may  be 
ubiquitous  within  the  surf  zone  and  are  believed  to  be 
related  to  forced  surf  zone  eddies. 

2.  Field  Site 

[7]  Observations  were  obtained  on  a  beach  with  quasi- 
periodically  spaced  rip  current  channels  at  Sand  City, 
Southern  Monterey  Bay,  California,  during  the  Rip  Current 
Field  Experiment  (RIPEX)-Steep  Beach  Experiment  (SBE) 
during  the  months  of  April  and  May  2001.  The  beach  face 
was  relatively  steep  with  a  slope  of  about  1:10  with  beach 
cusps  0(35  m),  which  flattens  out  to  a  low-tide  terrace 
(1:100)  with  quasiperiodic  O  (125  m)  incised  rip  channels, 
continuing  farther  offshore  with  a  1:20  slope  (Figure  2). 
Two  primary  arrays  were  deployed:  (1)  an  alongshore  array 
of  6  bidirectional  electromagnetic  (EM)  current  meters  that 
traversed  the  shore-connected  shoals  and  rip  channels,  and 
(2)  a  cross-shore  array  of  pressure  and  EM  current  meters 
on  the  shore-connected  shoal;  these  sensors  were  sampled 
continuously  and  synchronously  at  8  Hz.  EM  current  meters 
located  on  the  shore-connected  shoal  were  not  always 
submerged  at  lower  tides.  An  upward  looking  acoustic 
Doppler  current  meter  (ADCP)  was  deployed  offshore  of 
the  rip  channel  and  sampled  at  0.1  Hz. 

[8]  A  directional  wave  rider  buoy  located  ^650  m 
offshore  in  17  m  depth  measured  offshore  wave  conditions. 
Considerable  variation  in  wave  characteristics  occurred 
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Figure  2.  Bathymetry  for  year  day  117  obtained  by  combining  hydrographic  survey  system  mounted 
aboard  a  personal  watercraft  (high  tide)  and  a  backpack  mounted  to  a  person  walking  the  shoals  (low 
tide).  The  dots  indicate  pressure  sensors,  and  squares  represent  colocated  pressure  and  velocity. 


yearday 


Figure  3.  Wave  climate  measured  by  directional  wave  buoy  located  in  17  m  water  depth,  approximately 
650  m  offshore,  (top)  Hmo,  (middle)  peak  wave  period,  and  (bottom)  peak  wave  direction.  Shaded  regions 
represents  periods  of  measurement  discussion. 
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Figure  4.  Spectra  of  (top)  pressure,  (middle)  cross-shore,  and  (bottom)  alongshore  velocity  for  PUV 1 1 
located  within  the  rip  channel  for  year  day  124.1818.  The  regions  between  df  —0.004  Hz  have  a  DOF 
with  df  of  0.0001  Hz,  between  0.004  and  0.05  Hz  have  22  DOF  with  a  df  of  0.001  Hz  and  between  0.05 
and  0.1  Hz  have  42  DOF  with  a  df  of  0.0021  Hz.  The  region  to  the  left  of  the  vertical  line  at  0.004  Hz 
represents  VLFs. 
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time(yearday) 

Figure  5.  Time  series  of  (top)  pressure,  (middle)  cross-shore  velocity,  and  (bottom)  alongshore  velocity 
for  PUV  11  for  year  day  124.  The  black  line  is  the  raw  signal  and  the  white  line  in  the  low-pass  filtered 
signal  (FC  =  0.004  Hz)  representing  the  VLFs. 
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Figure  6.  Time  series  at  PUV11  for  (a,  b,  and  c)  year  day  111,  (d,  e,  and  f)  year  day  126,  and  (g,  h, 
and  i)  year  day  130.  The  black  line  is  the  raw  signal  for  pressure  (Figures  6a,  6d,  and  6g),  cross-shore 
velocity  (Figures  6b,  6e,  and  6h),  and  alongshore  velocity  (Figures  6c,  6f,  and  6i)  and  the  white  is  the 
low-pass  signal  (FC  =  0.004  Hz)  for  each  measurement.  Note  the  pressure  sensor  did  not  work  properly 
for  PUV 1 1  on  year  day  111. 


throughout  the  experiment  (year  days  100-140)  (Figure  3). 
In  spite  of  the  large  temporal  variability  in  wave  height  and 
period,  the  mean  wave  incidence  angle  was  predominantly 
shore-normal.  Beach  face  and  bathymetry  surveys  were 
obtained  using  a  differential  kinematic  global  positioning 
system  (KGPS)  mounted  on  a  sonar  equipped  personal 
watercraft  [ MacMahan ,  2001],  in  a  backpack  placed  on  a 
person,  and  on  an  all-terrain  vehicle. 

[9]  The  mean  flow  field  was  represented  by  cellular 
circulation,  with  predominant  shoreward  flow  on  the 
shore-connected  shoals  and  seaward  flow  through  the  rip 
channels  that  persisted  throughout  the  experiment,  except 
for  low  wave  conditions.  For  an  additional  discussion  of  the 
experiment,  see  MacMahan  et  al.  (submitted  manuscript, 
2004). 

3.  Very  Low  Frequency  Motions 

3.1.  Observations 

[10]  Similar  to  the  experience  encountered  by  Oltman- 
Shay  et  al  [1989],  we  often  experienced  significant  pulsa¬ 
tions  during  the  maintenance  and  inspection  of  the  EM 
current  meters.  Initially,  these  motions  were  attributed  to 


infragravity  rip  current  pulsations  [MacMahan  et  al .,  2004]. 
It  was  not  until  examining  the  time  series  of  the  current 
meter  data  that  energetic  very  low  frequency  motions 
(VLFs)  were  identified,  fluctuating  on  the  order  of  15  min. 
Owing  to  the  nearly  shore  normal  waves,  there  were 
minimal  (if  any)  longshore  currents.  Therefore  these 
motions  are  not  associated  with  the  shear  instabilities  of 
longshore  currents. 

[11]  To  obtain  high  resolution,  spectra  was  calculated 
from  2.28-hour  records,  resulting  in  frequency  resolution 
(df)  of  0.0001  Hz  with  2  degrees  of  freedom  (DOF).  The 
total  record  was  quadratically  detrended  to  remove  the  tidal 
signature.  A  significant  amount  of  energy  is  found  within 
the  VLF  band  for  both  the  cross-shore  and  alongshore 
velocity  components,  while  the  energy  within  the  sea 
surface  elevations  is  minimal  (Figure  4).  Note  that  energy 
is  significant  within  the  infragravity  (0.004-0.04  Hz)  and 
swell  and  sea  (0.04-0.35  Hz)  frequency  bands  for  sea- 
surface  elevation  and  cross-shore  velocities  and  minimal  for 
alongshore  velocities. 

[12]  Examples  of  VLFs  are  examined  in  the  low-passed 
(frequency  cutoff  (FC)  FC  =  0.004  Hz)  time  series  for 
pressure,  cross-shore,  and  alongshore  velocity  measure- 
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Figure  7.  Significant  wave  height  measured  at  the  wave  buoy  for  reference  is  shown  in  the  first  panel. 
Variance  of  UyLFrms  for  PUV6  located  within  the  rip  channel  is  shown  in  the  second  panel.  Tidal 
elevation  measured  at  the  NO  A  A  station  for  reference  is  shown  in  the  third  panel.  The  top  three  plots 
span  the  time  period  for  year  days  107-140,  while  the  bottom  three  plots  are  expanded  views  of  year 
days  120-130  to  examine  the  tidal  variability. 


ments  on  year  day  124  for  PUV11  located  within  the  rip 
channel  (Figure  5).  The  pressure  signal  has  minimal  VLF 
signal,  whereas  the  cross-shore  and  alongshore  velocities 
have  significant  VLF  signals  (Figure  5),  which  are  persis¬ 
tent  throughout  the  records. 

[13]  Different  VLF  behaviors  were  observed  during  the 
experiment.  Time  series  of  sea  surface  elevations,  cross¬ 
shore,  and  alongshore  measurements  for  sensor  PUV 1 1  for 
three  days  (year  days  111,  126,  and  130)  are  examined  in 
Figure  6.  For  year  day  111  the  measurements  illustrate  a 
segment  when  there  is  a  dramatic  increase  in  swell  energy 
and  the  VLFs  spin  up  (Figures  6a,  6b,  and  6c).  For  year 
day  126  the  alongshore  component  of  the  VLFs  are 
modulated  by  the  tide  with  increased  VLF  velocities 
occurring  at  lower  tides  (Figures  6d,  6e,  and  6f),  while 
year  day  130  demonstrates  a  period  when  the  VLFs  remain 
relatively  constant  throughout  the  day  (Figures  6g,  6h, 
and  6i).  For  all  three  days  the  alongshore  VLF  component 
qualitatively  appears  larger  than  the  associated  cross-shore 
component. 

3.2.  VLF  Statistics 

[14]  Low-passed  root-mean  squared  (rms)  VLF  values  for 
sea  surface  elevation,  cross-shore,  and  alongshore  velocity 
(hrms,VLF,  u rms, vlf,  H^vlf)  are  obtained  by  integrating  the 
energy  density  spectra,  E(f),  within  the  0-0.004  Hz  fre¬ 


quency  band.  Average  VLF  wave  period  (first  moment 
wave  period,  rm0i,vLF)  was  computed  instead  of  peak 
frequency  owing  to  the  sensitivity  of  the  measure  to  the 
frequency  resolution.  The  average  VLF  period  is  ~13  min 
for  both  cross-shore  and  alongshore  VLF  motions.  The 
distribution  of  VLF  periods  vary  slightly  between  locations. 
Most  VLF  periods  fall  within  the  500-1500  s  temporal 
range  with  mean  period  of  800  ±  300  s. 

[15]  The  temporal  and  spatial  velocity  variations  are 

examined.  {7rmSjVLF,  is  computed  from  +  'Wvi.f 

to  eliminate  any  potential  bias  associated  with  current 
meter  orientation.  t/nns,VLF  for  PUV6  has  daily  variability, 
but  the  overall  velocity  trend  increases  with  increasing 
sea-swell  energy  (Figure  7).  Qualitatively,  some  of  the 
daily  variability  can  be  associated  with  tidal  modulation, 
which  results  in  an  inverse  relationship  with  VLF  veloc¬ 
ities  increasing  during  decreasing  tidal  elevations.  This 
relationship  is  visibly  evident  for  year  days  126-130 
(Figure  7,  bottom  panel),  which  is  a  period  of  relatively 
constant  sea-swell  wave  energy.  VLF  rms  values  were 
computed  with  75%  overlap  to  provide  for  continuous 
measure  and  band-pass  filtered  to  retain  the  tidal  signa¬ 
ture,  resulting  in  an  overall  VLF  tidal  significant  (95%) 
correlation  of  r  =  0.28. 

[16]  Total  VLF  velocity  magnitudes  (t/rms,VLF)  were 
computed  for  year  days  111-120  to  examine  spatial  vari- 
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Figure  8.  f/rms,VLF  averaged  over  year  days  111-120  (top)  within  the  alongshore  array  and  (bottom)  in 
the  cross  shore.  The  95%  confidence  interval  is  represented  by  the  thickness  of  the  C/rms,VLF  curve.  The 
depth  profiles  are  shown  as  a  thick  curve  (depth  labeled  on  the  right-hand  side).  UrmsyLF  measured  at  by 
the  ADCP  in  the  neighboring  rip  channel  is  represented  by  black  circle.  The  dashed  line  represents  zero 
intensity  and  mean  sea  level. 
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Figure  9.  C/rms,VLF  plotted  versus  (left)  Ume an  measured  at  PUV6  and  (right)  Hmo. 
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Figure  10.  Cross  correlations  of  (top)  i/VLF  and  (middle)  vVLF  with  respect  to  measurements  at  PUV6 
located  within  the  rip  channel  at  y  =  0  in  the  alongshore  array  at  zero  time  lag.  Dashed  line  represents 
location  of  the  rip  channels,  (bottom)  Hmo  given  for  reference. 


ability,  which  is  representative  of  other  days  throughout  the 
experiment.  There  are  minimal  differences  in  the  alongshore 
transect  of  UrmsyLF  (Figure  8).  The  95%  confidence  interval 
is  the  thickness  of  the  line.  UrmsyLF  is  unaffected  by  the 
alongshore  bathymetric  variability  of  the  shore-connected 
shoals  and  rip  channels  (Figure  8)  and  is  not  focused  within 
the  rip  channels,  differing  from  observations  by  Haller  and 
Dalrymple  [2001].  This  suggests  that  the  observed  VLF 
motions  are  not  a  result  of  rip  current  instabilities  but  are 
developed  by  a  different  mechanism.  The  cross-shore  mag¬ 
nitudes  are  relatively  constant  within  the  surf  zone  but 
rapidly  decay  outside  the  surf  zone  (Figure  8).  Similar 
VLF  values  are  observed  for  measurements  outside  the  surf 
zone  of  the  shore-connected  shoal  and  the  rip  channel 
(Figure  8).  Thus  the  VLFs  are  energetic  (^0.2  m/s)  within 
the  surf  zone  with  minimal  alongshore  variability  and 
decrease  (^0.05  m/s)  outside  the  surf  zone. 

[17]  An  increase  in  rip  velocity  should  result  in  an 
increase  in  shear  as  the  rip  channel  width  is  constrained. 
If  the  VLF  motions  are  related  to  the  instabilities  of  the  rip 
currents,  then  an  increase  in  rip  current  velocity  shear 
should  increase  instability  intensities.  It  has  been  found  that 
rip  current  velocities  increase  with  increasing  wave  energy 
and  decreasing  tidal  elevation  [ Brander  and  Short ,  2001; 
Dronen  et  al .,  2002;  Haller  et  al .,  2002;  MacMahan  et  al., 
submitted  manuscript,  2004]  and  that  the  C/nns,vLF  qualita¬ 
tively  is  related  with  sea-swell  energy  and  tidal  elevation 
(Figure  7).  The  UrmsyFF  for  sensors  within  the  surf  zone  are 
plotted  correlated  with  the  mean  (90  min)  rip  current 


velocity  for  PUV6,  which  are  significantly  correlated  (r2  = 
0.18)  at  the  95%  confidence  level  (Figure  9).  Urrns?vLF  has 
a  greater  statistically  significant  (95%)  correlation  with  H0 
(r2  =  0.49)  than  with  Umean  (r2  =  0.18),  suggesting  that 
VLFs  are  influenced  by  external  wave  forcing. 

3.3.  VLF  Spatial  Behavior 

[is]  Cross-correlation  analysis  was  performed  to  deter¬ 
mine  alongshore  spatial  relationships  of  VLFs.  The  low- 
passed  (<0.004  Hz)  cross-shore  and  alongshore  velocities 
within  the  rip  channel  (PUV6)  are  correlated  with  adjacent 
sensors  in  the  alongshore  array  (PUV1,  PUV2,  PUV8, 
PUV9,  PUV1 1)  and  presented  at  zero  time  lag  as  time  stacks 
(Figure  10).  For  the  VLF  cross-shore  velocities  in  the 
alongshore  array,  the  correlation  varies  spatially  with  positive 
correlation  occurring  in  the  rip  channels  and  negative  corre¬ 
lation  occurring  on  the  shore-connected  shoals.  The  spatially 
varying  correlations  suggest  a  relationship  associated  with 
the  underlying  bathymetry  and  rip  current  flow  character¬ 
istics  of  vorticity  motion.  Similar  analysis  was  performed  on 
the  VLF  alongshore  velocities  (Figure  10,  middle  panel). 
High  positive  correlation  occurs  across  the  alongshore  array 
for  the  alongshore  velocities  during  energetic  wave  events. 
During  mild  wave  conditions  the  correlation  is  poor.  There 
was  ^800  s  time  lag  associated  with  the  correlation  maxima 
of  the  cross-shore  and  alongshore  VLF  velocities,  represen¬ 
tative  of  the  Tm0 i  vlf  of  the  measurements. 

[19]  Frequency  wave  number  ( f-ky )  spectra  were  calculated 
using  the  iterative  maximum  likelihood  estimator  (IMLE) 
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Figure  11.  F-ky  spectra  computed  using  IMLE  method  for 
(top)  cross-shore  and  (bottom)  alongshore  velocities  for  the 
alongshore  array  for  year  day  130.75.  The  dashed  lines 
represent  mode  0  and  1  edge  wave  dispersion  curves  based 
on  numerical  experiments. 

[ Pakwa ,  1983]  evaluated  by  Oltman-Shay  and  Guza  [1987] 
for  detection  of  edge  waves  and  later  by  Oltman-Shay  et  al. 
[1989]  for  the  detection  of  longshore  current  shear  insta¬ 
bilities.  A  basic  assumption  of  the  IMLE  is  spatial  homo¬ 
geneity,  which  is  not  necessarily  expected  over  bathymetry 
with  quasiperiodic  rip  channels.  The  IMLE  was  tested  and 
the  comparisons  of  IMLE  field  measurements  and  IMLE 
and  Fourier  model  simulations  suggest  the  IMLE  is  appli¬ 
cable  for  this  particular  beach  (Reniers  et  al.,  submitted 
manuscript,  2004)  through  comparisons  of  IMLE  field 
measurements  and  IMLE  and  Fourier  model  simulations. 
Edge  wave  dispersion  curves  for  this  beach  have  been 
determined  by  numerical  experiments,  see  Reniers  et  al. 
(submitted  manuscript,  2004). 

[20]  The  energy  within  the  f-ky  spectra  is  focused  within 
the  VLF  band  (<0.004  Hz)  and  nongravity  region  (Figure  1) 
for  both  the  alongshore  and  cross-shore  velocity  (Figure  11). 
This  is  unlike  the  characteristics  of  shear  instabilities  of 
longshore  currents,  where  energy  is  present  at  higher  fre¬ 
quencies  [ Oltman-Shay  et  al .,  1989].  For  the  cross-shore  f-ky 
spectra  (Figure  12,  middle  panel)  the  energy  within  the  VLF 


band  is  generally  concentrated  around  ky  =  ±0.008  m-1  (X^  = 
125  m),  approximately  equal  to  rip  channel  spacing,  while 
for  the  alongshore  f-ky  spectra  (Figure  12,  bottom  panel), 
energy  is  generally  focused  around  ky  =  0  m_1.  It  is 
surprising  that  there  are  two  different  spatial  scales  associ¬ 
ated  with  the  VLFs,  whereas  longshore  current  instabilities 
generally  have  similar  spatial  scales  for  both  cross-shore  and 
alongshore  f-ky  spectra  [ Oltman-Shay  et  al .,  1989]. 

[21]  F-ky  spectra  were  integrated  over  the  VLF  bands  for 
year  days  107-138,  and  plotted  as  a  wave  number  spectra 
time  stack  (Figure  12).  The  energy  within  the  VLF  cross¬ 
shore  velocities  is  generally  focused  around  rip  channel 
spacings,  ky  =  ±0.008  m_1,  and  higher  harmonics  of  the  rip 
channel  spacing  throughout  the  entire  experiment.  While 
the  energy  within  VLF  alongshore  velocity  is  focused  about 
ky=  0  m_1,  suggesting  a  wavelength  that  is  infinitely  long 
or  no  alongshore  propagation  (see  below).  The  energy 
increases  with  increasing  sea-swell  energy. 

4.  Discussion 

4.1.  Surf  Zone  Eddies  (SZEs) 

[22]  Peregrine  [1998]  conceptualized  the  behavior  of  free 
■|  eddies  within  the  surf  zone.  He  hypothesized  for  a  planar 
?  beach  that  a  single  eddy  would  travel  parallel  to  the  shore- 
^  line,  while  two  opposite-signed  eddies  (dipoles)  would 
~  induce  predominantly  cross-shore  motion.  He  also  suggested 

that  the  eddy  migration  is  limited  by  deeper  water  to 
approximately  two  surf  zone  widths.  When  more  than  two 
eddies  exist  in  the  surf  zone,  it  is  hypothesized  that  eddies 
will  propagate  randomly  in  the  cross  shore  at  various 
alongshore  locations. 

[23]  Ryrie  [1983]  showed  that  a  single  wave  group  packet 
approaching  the  shoreline  at  an  oblique  angle  with  a 
particular  alongshore  length  produces  surf  zone  eddies  (on 
the  order  of  surf  zone  scale),  which  are  slowly  dissipated  by 
bottom  friction.  A  wave  group  packet  has  a  particular  cross¬ 
shore  and  alongshore  spatial  scale.  Buhler  and  Jacobson 
[2001]  investigated  a  pair  of  oblique  wave  group  packets 
(dipole,  two  circulation  cells  of  opposite  signs),  which 
produced  eddies. 

[24]  Initial  observations  of  SZEs  were  documented 
during  the  nearshore  sediment  transport  study  by  Tang  and 
Dalrymple  [1989],  who  related  observations  of  migrating  rip 
currents  on  a  near  planar  beach  with  wave  group  forcing. 
This  observation  was  later  conceptualized  and  verified  in  the 
laboratory  by  Fowler  and  Dalrymple  [1990].  They  docu¬ 
mented  rip  currents  (SZEs)  propagating  slowly  at  speed,  C  zz 
Ak/Af,  forced  by  wave  groups  formed  by  two  different  wave 
trains  from  different  directions,  (Aky  =  Arising)  —  A;2sin(02) 
and  A/=/i  —  f2).  When  A ky  is  approximately  zero,  there  is 
no  alongshore  propagation  (Figures  11  and  12).  This  mech¬ 
anism  is  a  generalized  version  of  Dalrymple' s  [1975]  gen¬ 
eration  of  stationary  rip  currents  formed  by  two  intersecting 
(monochromatic)  wave  trains  at  equal  and  opposing  angles. 
As  Fowler  [1991]  increased  the  number  of  energetic  direc¬ 
tional  components,  the  hydrodynamic  response  in  the  f-ky 
spectra  became  broader. 

[25]  Reniers  et  al.  [2004]  generalized  these  previous 
approaches  (above)  with  realistic  directional  wave  spectra 
and  found  random  patches  of  wave  group  packets  entering 
the  surf  zone.  These  wave  group  packets  create  gradients  in 
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Figure  12.  Significant  wave  height  measured  at  (top)  PUV4  provided  for  reference.  Integrated  energy 
within  the  VLF  band  for  (middle)  cross-shore  and  (bottom)  alongshore  velocities  in  the  alongshore  array. 
Bold  lines  represent,  for  the  t/VLF  ky  time  stack,  the  approximate  rip  channel  spacing  (ky  =  ±0.008  m-1), 
while  for  the  UVLF  ky  time  stack,  ky  =  0  m-1. 


radiation  stresses  as  the  waves  break,  which  force  surf  zone 
eddies  that  have  preferential  scale  based  on  A ky  They  found 
a  morphodynamic  response  of  incised  rip  channels  at 
similar  scales. 

[26]  It  is  hypothesized  that  the  breaking  waves  of  group 
packets  induced  random  alongshore  gradients  in  radiation 
stresses,  forcing  SZEs  with  rip  currents  to  oscillate  within 
the  RIPEX  field  measurements.  However,  as  opposed  to  the 
SZEs  migrating  alongshore  within  the  surf  zone  [Fowler 
and  Dalrymple ,  1990]  or  migrating  freely  [ Peregrine , 
1998],  there  is  morphologic  feedback  occurring  that  couples 
the  rip  current  cells  to  the  underlying  bathymetry.  The 
random  wave  group  packets  force  the  cells  to  oscillate  within 
the  surf  zone  on  0(15  min)  timescales.  This  explains  the 
correlation  with  UrmsyLF  and  H0  as  well  as  the  large  amount 
of  scatter  (Figure  10).  The  complexity  of  this  beach  makes  it 
difficult  to  completely  evaluate  this  process.  The  SZE  rms 
magnitudes  are  spatially  homogenous  within  the  surf  zone 
(Figure  8).  Cross-correlation  analysis  with  the  cross-shore 
velocities  suggest  eddy  behavior  with  the  flows  being  out  of 
phase  between  the  rip  channels  and  shore-connected  shoals 
(Figure  10).  This  is  further  supported  by  the  f-ky  spectral 
estimates,  which  suggest  eddies  that  have  a  morphologic 
coupling  (Figures  11-12). 

[27]  The  spatial  scales  of  the  observed  f-ky  spectra  can  be 
explained  by  the  simple  schematic  in  Figure  13.  If  rip  current 


cells  exist  and  oscillate  within  the  VLF  band  in  the  cross 
shore,  this  will  produce  the  spatial  scales  observed  in  the  f-ky 
spectra  (Figures  11-12).  A  predominantly  cross-shore  dis¬ 
placement  (oscillation)  will  produce  energy  focused  around 
the  rip  current  spacings  for  cross-shore  velocities.  Taking  a 
slice  in  the  alongshore,  cross-shore  velocities  produce  a 


rip  current  rip  current 


Figure  13.  Hypothetical  rip  current  cells  at  t  =  0 
represented  by  solid  lines.  At  some  time  later  ( t  >  0),  rip 
current  cell  are  displaced  (dashed  line).  This  motion  induces 
alongshore  spreading  around  ky  =  0  m-1  for  the  alongshore 
velocities  in  the  f-ky  spectra,  while  producing  energy  at  the 
rip  channel  spacing  for  ky  =  ±0.008  m-1  for  the  cross-shore 
velocities. 
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Figure  14.  (top)  The  mean  flow  structure  of  the  conceptual  model  of  rip  current  velocity  cells  based  on 
equations  (l)-(3).  (bottom)  The  variance  of  the  velocity  field.  Solid  line  is  the  location  for  computations 
of  the  two-dimensional  FFT  representative  of  the  location  of  the  alongshore  array  in  the  field. 


sinusoidal  curve  which  oscillates  at  a  particular  frequency 
(VLF),  and  it  is  relatively  fixed  in  space,  which  is  associated 
with  the  rip  current  spacings.  The  alongshore  velocities  will 
remain  relatively  the  same,  except  there  will  be  a  slight 
alongshore  displacement,  which  smears  the  energy  around 
ky  =  0  m-1  (Figures  11-12).  This  indicates  that  there  is 
morphological  feedback  coupled  with  the  rip  current  eddy 
cells,  as  they  are  topographically  controlled  by  the  shore- 
connected  shoals  and  rip  channels,  but  there  is  relaxation 
within  the  hydrodynamics  to  allow  the  rip  current  cells  to 
move.  Thus  the  SZEs  at  RIPEX  are  being  forced  by  radiation 
stress  gradients  associated  with  breaking  of  wave  groups  and 
are  not  free  eddies  traveling  along  the  coast. 

4.2.  Kinematic  Conceptual  Model 

[28]  A  conceptual  model  of  rip  current  cell  circulations  is 
developed  to  aid  in  interpreting  the  SZE  behavior  observed 
in  the  f-ky  spectral  observations  (see  Figures  11  and  12).  In 
the  model  a  two-dimensional  Gaussian  stream  function  is 
utilized  to  describe  the  cells,  defined  by 

V(x,y,  t)  =  yy  1  (!) 

n=  1 

where  N  is  the  number  of  cells,  —\n  makes  the  field 
spatially  periodic,  A  is  a  constant  amplitude  of  the  Gaussian 
field,  X  is  the  length  of  the  cell,  xc,  yc  represents  the  center 
of  the  cells,  and  t  is  time.  Assuming  a  rigid  lid,  velocities  are 
defined  as  stream  function  derivatives, 


Each  rip  current  cell  is  converted  into  velocity  vectors. 
Multiple  rip  current  cells  are  located  in  space,  representing 
the  approximate  locations  of  rip  currents  in  the  RIPEX 
bathymetry  (Figure  2).  The  overlapping  of  neighboring 
rip  currents  cells  is  additive,  such  that  the  section  of  cells 
located  in  close  proximity  are  faster  than  cells  with  sections 
which  are  further  apart  (superposition)  (Figure  14  (top)).  Rip 
current  cells  are  separated  in  the  vicinity  of  onshore  transport 
(shore-connected  shoals)  and  close  together  for  offshore 
transport  (rip  channels),  producing  faster  offshore  directed 
flows  (representative  of  rip  currents).  The  model  has  regions 
of  zero  velocity,  which  would  not  be  the  case  in  nature. 

[29]  The  center  (xc,  yc )  of  each  rip  current  cell  is  sinusoi¬ 
dally  oscillated  in  time, 

12 

xc{n,i)\yc(n,t)\  =  ^dXk[dYk\  sin^Ttf  +  e*),  (3) 

k= 1 

with  uniform  random  distributions  of  amplitudes  (dXh  dY k), 
frequencies  (0  <  fk  <  0.004  Hz),  and  phases  (e*),  which 
cause  the  rip  current  cells  to  oscillate  in  the  cross-shore  and 
alongshore  directions.  Here  dX  and  dY  are  the  potential 
maximum  amplitudes  of  the  uniform  distributions.  Simula¬ 
tions  were  run  for  approximately  two  hours  with  a  sampling 
interval  of  5  s,  similar  to  field  velocity  measurements  for  the 
f-ky  spectra.  F-ky  spectral  estimates  are  computed  using  two- 
dimensional  (2-D)  fast-Fourier  tansforms  (FFTs)  with  a 
hamming  window  resulting  in  40  DOF  with  df=  0.0004  Hz 
for  measurements  along  y  =  88  in  Figure  14  (solid  line). 
F  —  ky  spectra  represent  the  fluctuations  associated  with  the 
system,  not  the  mean  flow  structure  of  the  rip  current  cells 
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Figure  15.  F  —  ky  spectral  estimates  using  2-D  FFT  of  kinematic  conceptual  models  of  rip  current 
cell  oscillations;  (first  row)  oscillations  of  cross-shore  rip  current  cell  dX  =  15  m,  dY  =  0  m,  (second 
row)  oscillations  of  alongshore  rip  current  cell  dX  =  0  m,  dY  =  1 5  m,  (third  row)  cross-shore  and 
alongshore  oscillation  of  rip  current  cells  dX=  15  m,  dY=  4  m,  (fourth  row)  IMLE  for  the  cross-shore  and 
alongshore  oscillation  of  rip  current  cells  dX=  15  m,  dY=  4  m,  (fifth  row)  2-D  FFT  for  VLF  amplitude, 
A(t),  modulations  of  rip  current  cells. 


(Figure  14).  Three  scenarios  were  numerically  performed: 
(1)  rip  current  cells  that  only  oscillate  randomly  in  the  cross¬ 
shore  (dX=  15  m,  dY=  0  m),  (2)  rip  current  cells  that  only 
oscillate  randomly  in  the  alongshore  (dX=  0  m,  dY=  15  m), 
and  (3)  rip  currents  cells  that  oscillate  randomly  in  both 
directions  (dX  =  1 5  m,  dY  =  4  m). 

[30]  The  results  indicate  that  when  the  rip  current  cells 
only  oscillate  in  the  cross  shore,  the  ky  of  the  cross-shore 
velocities  is  equivalent  to  the  rip  current  spacings,  while 
the  ky  of  the  alongshore  velocities  is  spread  around  zero 
(Figure  15  (top)).  This  is  similar  to  the  measurements 
obtained  from  the  alongshore  array  of  the  field  experiment 
(Figure  11).  For  strictly  alongshore  oscillations  (Figure  15 
(second  panel)),  ky  is  focused  around  rip  current  spacing 
for  both  cross-shore  and  alongshore  velocities.  A  scenario 
with  equal  cross-shore  and  alongshore  movements  was 
also  performed  (not  shown),  which  produced  results  sim¬ 
ilar  to  the  case  with  just  alongshore  oscillations  (Figure  15 
(second  panel)).  The  last  case  is  for  a  slight  amount  of 


alongshore  oscillation  (4  m)  dominated  by  cross-shore 
oscillation  (15  m)  and  tends  to  smears  the  f-ky  spectra 
(Figure  15  (third  panel)),  which  are  most  comparable  to 
the  field  measurements  (Figure  11).  The  IMLE  method 
was  performed  on  the  last  case  with  a  spatial  array  similar 
to  the  field  array  (Figure  1)  with  some  additional  locations 
that  span  the  entire  alongshore  domain  (Figure  15  (fourth 
panel)).  The  IMLE  has  similar  results  to  that  of  the  2-D 
FFT  but  with  more  localized  energy  within  ky.  This 
suggests  a  plausible  description  that  during  the  experiment 
the  rip  current  cells  randomly  oscillated  in  the  cross  shore 
with  smaller  alongshore  oscillations. 

4.3.  SZE  Wave  Group  Behavior 

[31]  Another  potential  mechanism  that  can  explain  the 
SZE  motions  is  that  slow  modulations  in  wave  group  energy 
increase  the  velocities  of  the  rip  current  cells  through 
pressure  and  radiation  stress  gradients  [Bowen,  1969].  The 
randomness  of  the  wave  groups  produce  durations  of  large 
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Figure  16.  Contours  (color  scale  at  the  right)  of  cross-correlation  functions  for:  (first  row)  VLF  offshore 
wave  group  energy  (PUV4)  and  rip  channel  cross-shore  VLF  velocity  (PUV11);  (second  row)  VLF 
offshore  wave  group  energy  (PUV4)  and  rip  channel  alongshore  VLF  velocity  (PUV11)  as  a  function  of 
lag  and  time;  (third  row)  Hmo  measured  at  the  buoy  versus  time;  and  (fourth  row)  tidal  elevation  versus 
time. 


energy  followed  by  durations  with  less  energy,  which  would 
be  expected  to  induce  a  slowly  varying  rip  current  velocity 
field.  This  would  appear  as  though  the  rip  current  cells  were 
oscillating  in  the  cross  shore  as  described  by  the  conceptual 
model.  To  examine  this  mechanism,  a  modulation  of  Gauss¬ 
ian  amplitude  (A)  is  applied  to  equation  (1),  keeping  the 
center  of  the  cells  fixed,  where 

A{t)  =  A0  +  A'  sin(W  +  e),  (4) 

where  A0  is  a  constant,  A'  is  the  modulation  amplitude,  u;  is 
the  radian  frequency  with  the  VLF  band,  and  e  is  the  phase 
of  the  system.  The  VLF  amplitude  modulation  results  in 
energy  focused  around  the  rip  channel  spacings  for  both 
cross-shore  and  alongshore  velocities  in  f-ky  space  (Figure  15 
(bottom  panel)).  Though  the  energy  is  larger  in  the  cross¬ 
shore  velocity,  the  spatial  pattern  does  not  match  the  field 
observations. 

[32]  Cross-correlation  functions  were  computed  between 
the  VLF  wave  group  energy,  E(t),  at  PUV4  and  UVLF 
within  the  rip  channel  at  PUV11,  plotted  as  a  time  stack 
(Figure  16).  There  is  a  positive  correlation  occurring  at 
~50  s,  which  is  theoretical  time  for  wave  group  energy  to 
propagate  from  PUV4  to  PUV11  [MacMahan  et  al ., 
2004].  VLF  modulations  of  wave  group  energy  are  corre¬ 
lated  for  the  slow  oscillations  of  the  rip  current  cells. 
However,  there  is  no  clear  trends  associated  with  the 
correlation,  in  particular  during  times  of  increased  wave 
energy,  when  the  SZE  energy  increases.  The  temporal 
variation  of  wave  group  energy  is  only  part  of  the 
mechanism;  the  influences  of  the  A ky  also  induce  random 
eddies  into  the  surf  zone  at  VLF  frequencies  [Reniers  et 
al .,  2004].  Thus  it  is  the  spatial  and  temporal  aspects  of 
the  wave  groups  and  their  coupling  with  the  underlying 


bathymetry,  which  force  the  SZE  oscillations,  as  explained 
by  the  kinematic  conceptual  model. 

4.4.  Jet  Instability 

[33]  In  addition  to  the  rip  current  cell  oscillations  (de¬ 
scribed  above),  jet  instabilities  may  be  associated  with 
smaller  surf  zone  eddies.  In  a  laboratory  experiment,  Haller 
and  Dalrymple  [2001]  found  that  topographically  controlled 
rip  currents  [Haller  et  al ,  2002]  forced  by  monochromatic 
waves  were  often  unstable.  They  obtained  good  compar¬ 
isons  of  their  observations  with  a  model  based  on  the  linear 
instability  theory  of  a  narrow  jet.  Chen  et  al.  [1999] 
obtained  good  agreement  with  the  same  laboratory  measure¬ 
ments  of  rip  currents  using  the  fully  nonlinear  Boussinesq 
formulation.  Yu  and  Slinn’s  [2003]  numerical  evaluation  of 
wave-current  interactions  of  rip  current  flows  indicates  that 
the  rip  current  instabilities  originate  within  the  feeder 
channels  and  migrate  offshore. 

[34]  Some  of  the  laboratory  observations  by  Haller  and 
Dalrymple  [2001]  resemble  the  field  observations.  How¬ 
ever,  in  the  laboratory,  the  instabilities  are  concentrated 
within  the  rip  channel  neck  and  decrease  in  intensity  over 
the  neighboring  bars,  while  in  the  field  the  SZEs  are  of 
similar  magnitude,  regardless  of  sensor  location  within  the 
surf  zone  (Figure  8).  This  may  be  a  result  of  the  relative 
close  proximity  of  rip  current  cells  in  the  field  ( \r/wr  =  2) 
(Figure  2),  which  were  not  completely  independent,  whereas 
the  rip  current  systems  in  the  laboratory  were  relatively 
independent  (\r/wr  =  5),  where  \r  is  the  rip  channel  spacing 
and  wr  is  the  rip  channel  width.  The  characteristics  of  the 
morphology  are  also  different.  The  laboratory  experiment 
morphology  is  a  longshore-trough  bar  with  straight  rip 
channels  cutting  through  the  bar,  which  may  confine  the  flow 
more  than  the  complex  morphology  of  a  shore-connected 
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shoal  with  smaller  feeder  channels  and  incised  rip  channels 
that  increase  in  width  as  a  function  of  cross-shore  distance. 

[35]  The  Haller  et  al.  [2002]  laboratory  measurements  are 
scaled  using  an  undistorted  Froude  number  (Fr  =  0.2 -0.6) 
with  a  length  scale  ratio  of  1/50,  which  corresponds  to  a  rip 
current  system  in  nature  having  450  m  rip  spacing,  90  m  rip 
width,  breaking  wave  heights  1.3 -3. 8  m,  wave  periods  of 
5.7-7. 1  s,  and  mean  rip  current  velocities  of  1.0 -2.1  m/s. 
This  is  different  than  the  field  conditions  during  the  RIPEX 
field  experiment.  During  RIPEX  the  rip  channel  system 
corresponded  to  a  rip  current  spacing  of  125  m,  60  m  rip 
current  width,  wave  heights  at  P5  of  0.3-2  m,  4-10  s 
wave  period,  0-1  m/s  mean  rip  current  velocities,  resulting 
in  Fr  =  0-0.45. 

[36]  The  largest  difference  between  the  field  and  labora¬ 
tory  scaled  rip  currents  is  the  flow  speeds,  which  may  be  a 
result  of  the  differences  in  geometry,  rip  channel  alongshore 
spacing,  width  of  channel,  width  of  bar,  and  depth  of 
channel.  This  results  in  significant  differences  in  rip  current 
shear,  computed  as  the  mean  velocity  within  the  rip  channel 
divided  by  the  width  of  the  rip  channel.  The  shear  for 
the  laboratory  was  0.11  m/s/m  for  a  0.2  m/s  rip  flow  and  a 
1.82  m  rip  channel  width,  while  in  the  field  during  RIPEX 
the  shear  was  0.008  m/s/m  for  a  0.5  m/s  rip  flow  and  a  60  m 
rip  channel  width.  This  order  of  magnitude  difference 
suggests  that  the  mechanisms  responsible  for  the  instabil¬ 
ities  may  be  different  as  well  as  their  spatial  occurrences. 
Furthermore,  the  wave  forcing  is  also  different,  as  the 
laboratory  experiment  used  monochromatic  waves  and  field 
conditions  consisted  of  directionally  broad  wave  groups. 
Thus  the  field  observations  of  VLF  SZEs  do  not  appear  to 
be  associated  with  jet  instabilities. 

5.  Summary  and  Conclusion 

[37]  Energetic  motions  at  the  very  low  frequency 
(<0.004  Hz)  were  observed  throughout  the  RIPEX  experi¬ 
ment.  SZEs  were  present  everywhere  on  the  beach  system, 
composed  of  shore-connected  shoals  with  incised  rip 
channels.  SZE  velocities  were  ~25  cm/s  within  the  surf  zone 
with  alongshore  homogeneity  and  decreases  to  ^5  cm/s 
outside  of  the  surf  zone.  The  SZE  cross-shore  velocity 
energy  is  concentrated  within  the  very  low  frequency 
band  (<0.004  Hz)  and  with  alongshore  wave  number 
corresponding  to  the  morphodynamic  scale  of  rip  channels, 
differing  from  observations  of  the  energy  of  shear 
instabilities  of  the  longshore  current  residing  at  higher 
frequencies  [ Oltman-Shay  et  al ,  1989].  Insignificant  mean 
longshore  currents  were  present  during  the  field  measure¬ 
ments  to  generate  longshore  current  instabilities.  Previous 
field  measurements  of  rip  currents  seem  to  have  documented 
SZEs,  but  a  conclusive  interpretation  was  not  available 
[Smith  and Largier ,  1995;  Brander  and  Short ,  2001].  Haller 
and  Dalrymple  [2001]  illustrated  through  detailed  measure¬ 
ments  in  a  laboratory  and  an  analytical  model  that  the 
jet-like  flow  of  the  rip  current  can  be  unstable.  These  field 
measurements  of  SZEs  differ  from  the  laboratory  measure¬ 
ments  by  morphological  and  hydrodynamic  scales,  \r/wr, 
Fr ,  and  a  order  of  magnitude  difference  in  shear. 

[38]  SZE  energy  increases  with  increasing  sea-swell 
energy  and  with  some  possible  tidal  modulation.  There  is 
a  weaker  statistically  significant  correlation  between  SZE 


velocity  variation  and  Ume an  (Figure  9).  A  better  correlation 
exists  between  Urms?VLF  and  H0  (Figure  9),  indicating  that 
wave  forcing  may  be  responsible  for  triggering  the  rip 
current  cell  oscillation.  The  alongshore  velocities  are  cor¬ 
related  in  the  alongshore  when  the  SZEs  are  most  energetic 
(Figure  10),  while  the  cross-shore  velocities  are  out-of¬ 
phase  between  the  rip  channels  and  the  shore-connected 
shoals  (Figure  10).  F-ky  spectra  indicate  that  the  SZEs  of 
the  cross-shore  velocities  are  coupled  to  the  spacing  of 
the  rip  current  channels  (ky  =  ±0.008  m-1),  while  along¬ 
shore  velocities  are  concentrated  around  ky  =  0  m-1 
(Figures  11-12).  A  simple  kinematic  conceptual  model 
(Figure  13)  of  the  entire  rip  current  system  oscillating 
produced  results  similar  to  the  field  observations 
(Figure  14).  The  best  comparison  were  formed  when  the 
rip  current  cells  tend  to  oscillate  more  in  the  cross  shore 
than  in  the  alongshore  (Figures  11,  12,  and  15). 
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